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Q Abstract: We investigate the behavior of the retarded Green's function of a U(l) gauge 

r^i field in holographic A^ = 4 Super Yang-Mills plasma, taking the leading strong coupling 
corrections into account. First, we use the thermal limit of this quantity to determine the 

►5^ flow of the photon quasinormal mode spectrum away from the infinite 't Hooft coupling 

■^ limit, and after this specialize to a specific model of holographic thermalization, in which 

.^ we evaluate the corresponding spectral density. In the latter case, our primary interest lies 

CO in the pattern, with which the spectral density approaches its equilibrium form, and how 

■^ this process depends on the value of the coupling as well as the photon virtuality. All of the 

#y^ results obtained point consistently towards the weakening of the usual top/down pattern 

^~~i of holographic thermalization, once the coupling is decreased from the A = 00 limit. 
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1 Introduction 

Quantitatively understanding the processes driving the complicated field dynamics during 
a relativistic heavy ion collision presents an exceptionally difficult and topical challenge to 
QCD theorists. According to standard lore, the early onset of hydrodynamical behavior as 
well as the apparently very low value of the shear viscosity of the produced matter, both 
inferred from RHIC and LHC data, point towards the system being a strongly coupled, 
nearly ideal fluid [1, 2]. This has led to gauge/gravity methods becoming standard tools 
in the description of quark gluon plasma (QGP) physics, filling in gaps left by standard 
approaches such as perturbation theory and lattice simulations (see e.g. [3]). The success 
of the holographic approach relies largely on the fact that even though in their respective 
vacuum states QCD and conformal N = A Super Yang-Mills (SYM) theory are clearly very 
different, at high temperatures they both describe deconfined plasmas with qualitatively 
similar properties. 

A particularly useful recent development in the holographic description of QCD matter 
has been the mapping of the off-equilibrium dynamics driving thermalization in a strongly 
interacting plasma to the gravitational problem of black hole formation in asymptotically 



Anti de Sitter (AdS) spacetimes. Although this dual system is computationally almost 
equally challenging to approach, gauge/gravity calculations have already led to impressive 
results, such as a qualitative resolution of the famous fast apparent thermalization puzzle 
by demonstrating that the applicability of hydrodynamics does not require the system to be 
isotropic, let alone thermal [4-6]. Other recent applications of holography in this context 
include e.g. progress in the evaluation of off-equilibrium Green's functions [7] , developments 
in various thermalization scenarios (see e.g. [8-11]), as well as the determination of a host 
of kinetic theory parameters at strong coupling [12-15]. 

A particularly important challenge in the study of the thermalization process in heavy 
ion collisions is to identify the pattern, with which plasma constituents with different en- 
ergies approach their final thermal distributions. On the weak coupling side, a classical 
result identified elastic scatterings as the dominant process behind a bottom-up type ther- 
malization [16], which has recently received confirmation from classical gauge theory simu- 
lations [17] (see also [18] for a competing/complementing proposition emphasizing the role 
of plasma instabilities). In contrast, all holographic studies of thermalization in infinitely 
strongly coupled gauge theory have consistently pointed towards a top-down pattern, in 
which the modes with highest energy are the first to thermalize. This observation, which 
can be understood via a simple causal argument, indicates the probable existence of a 
transition between the two thermalization patterns at intermediate coupling. While the 
details of this conclusion certainly depend on the initial conditions the system is subjected 
to, it clearly poses an interesting topic for further study. 

A convenient way of studying the thermalization pattern of a plasma with given initial 
conditions and coupling strength is through the determination of various off-equilibrium 
Green's functions, whose evolution gives direct information about the rate with which 
plasma constituents of different energy approach the thermal limit. Due to technical com- 
plications inherent in their holographic determination (see e.g. [7]), one is in practice forced 
to work within the simplest and most symmetric models of thermalization, such as one 
where the process is dual to the radial collapse of a thin spherical shell in AdSs space and 
the subsequent formation of a black hole. This model, originally proposed by [19, 20] and 
further studied e.g. in [21], corresponds to an initial state rather far away from the early 
stages of a real life heavy ion collision. At the same time, it however allows for rather 
complicated calculations to be carried out, such as the determination of off-equilibrium 
spectral densities and particle production rates [21-24]. 

In the phenomenological description of the thermalizing plasma, one particularly in- 
teresting observable is the spectrum of emitted photons, as these particles interact with the 
plasma constituents so weakly that they escape the system almost unaltered. In thermal 
equilibrium, their production has been considered to leading as well as next-to-leading or- 
der (LO and NLO) in perturbative QCD [25, 26], to LO in the limits of weakly and strongly 
coupled AA = 4 SYM theory [27] , and in the latter theory even to NLO in a strong 't Hooft 
coupling expansion [28, 29], including corrections of order 1/A'^". The LO holographic 
calculations have since then been generalized to an out-of-equilibrium setup [22, 23], and 
later to NLO for the case of real photons [30]. These studies revealed interesting details 
of the virtuality and coupling strength dependence of the thermalization process, some 



of which were in fact rather striking. In particular, the NLO study of photon emission 
in [30] found indications of the top-down pattern of holographic thermalization shifting 
towards bottom-up type behavior, once the leading finite coupling corrections were taken 
into account. 

In the paper at hand, our aim is to continue the study of holographic thermalization 
through a further analysis of photon production in AA = 4 SYM plasma. In particular, 
we plan to fill a number of gaps left by our earlier works [22, 23, 30] by completing the 
following three tasks: 

1. Determine the quasinormal mode (QNM) spectrum of the photons from the corre- 
sponding retarded Green's function in thermal equilibrium, working at NLO in a 
strong coupling expansion. Study the flow of the QNM poles with A for several 
different virtualities. 

2. Compute the photon spectral density for different virtualities and couplings, limited 
only by the convergence of the strong coupling expansion. 

3. Perform an analytic WKB-type calculation for the asymptotic large-w behavior of 
the photon spectral density, accounting for the 0{l/\'^''^) corrections. 

We anticipate that these three results will shed more light on the conclusions of [30] , either 
strengthening or challenging the claims made therein. Doing this, we also hope to identify 
the most fruitful directions of future research within these topics. 

Our paper is organized as follows. Section 2 contains a brief introduction to the 
collapsing shell model of [19], while section 3 lays out the general strategy, with which 
we determine the photon Green's functions. In both of these cases, we begin from the 
A = oo limit, and in separate subsections consider the modifications introduced by including 
the leading 0(1/A^'^) corrections in the calculations. After this, section 4 is devoted 
to presenting the results for both the QNM spectra and the spectral densities, while in 
section 5 we draw our conclusions. Finally, the computational details related to the WKB 
calculation, discussed also in section 4, are left to appendix A. 

2 The setup 

2.1 Generalities 

We wish to use the collapsing shell model of [19, 20] to gain insights into the thermalization 
process of strongly coupled AA = 4 SYM plasma. On the dual gravity side, this model 
describes the gravitational collapse of a spherically symmetric, infinitely thin shell of matter 
in the radial direction of AdS5 space — a process dependent on the initial condition of the 
shell as well as its equation of state. Leaving these details unspecified, we parameterize the 
progress of the system by the radial location of the shell r^, which eventually approaches 
a value corresponding to the Schwarzschild horizon of the forming black hole, r/j. On the 
field theory side, the setup may be thought of as resulting from an instantaneous, spatially 
homogeneous injection of energy into the vacuum state of the theory; for details of the 
realization of such a system, see e.g. [31-33] 



Throughout our calculations, we work in the quasistatic limit of the dynamics, in 
which the motion of the shell is considered slow in comparison with the other time scales 
of interest. As discussed in [22], in the case of a shell composed of null dust (obeying the 
relation ds^ = 0), its approach towards the horizon takes the form 

= e -h , (2.1) 



where t denotes the time coordinate on the boundary (and thus in the dual field theory) 
and Tfi = l/(47rT), with T denoting the temperature of the field theory system in its final 
thermal state. The approximation of neglecting the inverse time scale associated with the 
motion of the shell in comparison with the frequencies (momenta) of the Green's functions 
of interest thus improves with time, and only breaks down in the limits of early times and 
very small energies. 

Following Birkhoff 's theorem, we know that outside the shell the metric of our system 
corresponds to that of an AdS black hole, whereas inside the shell we have pure AdSs 
space. Changing the radial coordinate to n = r^/r^, this means the metric obtains the 
compact form 

2 r 2 

ds"^ = -^ (f{u)dt'^ + dx"^ + dy^ + dz^) + ^ „^, . du^ , (2.2) 



where 




'-•> ^ 1 tor u > u^ . 



Here and in the following, we have chosen the subscript '-' to denote the inside and '+' 
the outside of the shell. As discussed in [22], the discontinuity of the time coordinate at 
the shell implies that frequencies measured inside and outside the shell must also differ, 
and are in fact related through 

iO- = ^^^, /„ = U{us). (2.4) 

Note that consistently with this notation, uj will below always refer to a frequency measured 
outside the shell. 

Finally, we note that a further issue resulting from the presence of the shell is the 
changing boundary conditions of various fields at r = r^ or u = Us. For an electric field E, 
a straightforward application of the Israel junction condition gives (for more details, see 
e.g. [21, 22]) 



(2.5) 

duE^{u,U;^)\u, = fmduE-i-{u,UJ+)\us , 

which will be an important ingredient of our determination of the retarded Green's function 
of a photon in the following sections. 



2.2 Finite coupling corrections 

In order to perform calculations beyond the usual A = oo limit of the AdS/CFT corre- 
spondence, we must account for all corrections of order a'^ in the inverse string tension. 
In particular, this requires supplementing the usual type IIB supergravity action by terms 
proportional to 7 = |C(3)A~2 ^ 

5iib = 50b + 75i\b + O(7'), (2.6) 

where the expansion is performed around the 0(7'') action 






Rio-\{dct>f-^{F,f. 



(2.7) 



and the correction term obtains the form [34, 35] 



L^ /■ ,10 ^ ^ ./^ X4 



In these expressions, Rio denotes the Ricci scalar, (p the dilaton field and F5 the five- form 
field strength, while C stands for the Weyl tensor and the tensor T is given by 

Ti. ., ^ — iV F^ -\- — ( F~^ /?+»""_ ■?/?+ /?+ »""^ (2 Q) 

labcdef — ''^ a^^cdef ^ ^q \^ abcmn^ def '^^ abfmn^ dec ) ' y"^'^) 

In the last equation, the index triplets {a, b, c} and {d, e, /} are understood to be first 
antisymmetrized with respect to all permutations, and the two triplets then symmetrized 
with respect to the interchange abc o def. One should also note that our notation for the 
various contractions in eq. (2.8) is only schematic, and e.g. the C^ term in fact denotes the 
combination 

'-^hmnk^pmnq^h ^ rsk' ^ '^pqmn'^h '^rsk' K"^-^^) 

For further details of this construction, see e.g. [28]. 

The derivation of the 7-corrected background metric from the above action is in prin- 
ciple a rather laborious task, which is however significantly simplified by the observation 
that the tensor T vanishes when evaluated with the 0(7*^) metric of eq. (2.2) [35]. This 
implies that the C^ term of eq. (2.10) is the only one giving nonzero contributions to the 
0(7) terms in the metric, eventually leading us to the result [34, 36, 37] 

ds^ = ^ i-f+iu) K^(u) dt^ + dx^) + , ^l , , P^(u) du^ + L'^iu) d^l , (2.11) 
u 4u^j_l_(u) 

where the various functions read 

K{u) = eT [«(")+4b(«)] ^ p(^) = e7b(«) ^ _^(y) = e7c{«) ^ 
1625 2 .^r 4 10005 6 
8 lb 

325 2 1075 4 4835 . 

h{u) = u^ -\ u* u" , 

^ ^ 8 32 32 ' 

1 5 

c{u) = -{l + u^)u\ (2.12) 

The 7-corrected relation between r^ and the field theory temperature finally obtains the 
formr;,=7rr/(l+26|7). 



3 Photon production 

3.1 Generalities 

Next, we wish to apply the falhng sheh setup to an investigation of photon production 
during the thermahzation process, complementing the treatment of [22, 23] by keeping the 
virtuality of the photons v = {uj^ — q^^juP' arbitrary. To this end, we denote the photon 
four- momentum by Q = (w, q) and parameterized the magnitude of the three- momentum 
q by g = cw, letting c = \/\ — v vary from to 1. These two extreme cases correspond 
respectively to dileptons at rest, considered in [22], and to real photons, cf. [23]. Again, we 
begin our treatment from the A = oo limit, continuing only in the next subsection to the 
effects of 0(7) corrections. 

In the context of AA = 4 SYM, by 'photon' we mean a gauge field coupled to a conserved 
current corresponding to one of the U(l) subgroups of the SU(4) R-symmetry of the theory, 
chosen so that all the AA = 4 scalars and fermions obtain equal charges. The quantity we 
are after is then the retarded Green's function of this field, denoted by Il^iy{Q), from the 
transverse and longitudinal components of which one obtains the trace of the corresponding 
spectral density, 

X^^iQ) = -UmU^iQ) - 2Imn||(Q) . (3.1) 

Using the fact that the fluctuation dissipation theorem is known to hold in the quasistatic 
limit [22], this allows us to evaluate the photon Wightman function 

U<,{Q) = n{u)Xf.u{Q) , n{oo) = l/{e^^ - 1) , (3.2) 

which is further related to the differential photon production rate 

dT ^ «<'^n<,(Q) 

with a denoting the fine structure constant. 

The retarded Green's function is a convenient quantity for holographic calculations, 
as it can be evaluated with the (by now rather standard) methods developed in [38]. To 
this end, we begin from the classical equation of motion (EoM) of an U(l) vector field in 
curved spacetime, 

^ =5^[\/=55^^5""i^pa] = , Fp„ = dpA^ - d^Ap , (3.4) 



y-9 

from which we obtain differential equations for the transverse and longitudinal components 
of the corresponding electric field. In the metric of eq. (2.2), this produces the momentum 
space EoMs (see also [27, 39]) 



/'(n) u^ - ff ju) 

--— S_L U) + — - 

/(u) uf{u) 






where f{u) denotes the function defined in eq. (2.3) and where we have introduced the 
dimensionless parameters u = w/{2TrT), q = g/(27rT). Already from here, two special 
limits are clearly visible: For deeply virtual photons at q = (dileptons at rest) the two 
EoMs agree, while for real photons with q = Co only the transverse component contributes. 
In these two cases, it thus suffices to consider only eq. (3.5). 

Next, we proceed to solve the above equations, taking as an example the transverse 
electric field. Outside the shell, u < Us, we have f{u) = /+(«), causing the equation (3.5) 
to have a singular points at the horizon, n = 1, with the indicial equation 



q" + 



UJ 



0. 



(3.7) 



The two solutions of this equation, a = ^iu/2, are readily identified as the infalling and 
outgoing modes, of which the former, corresponding to the minus sign, is the physical one in 
the thermal limit. In the presence of a falling shell, we must however keep both functions, 
and thus write the general 'outside' solution as a linear combination of the infalling and 
outgoing modes 

E+{u) = c+Ei4u) + c-Eout{u) , (3.8) 

with c± denoting unknown coefficients. The two functions appearing on the right hand 
side are finally solved numerically using the ansatze 



E^ 



i/out{u) = (1 - U)^ 2 yin/out{u) 



(3.9) 



with the functions yin/out('") satisfying the boundary conditions yin/out(^ = 1) = 1 and 

y:n/out(^=l)=0- 

Having obtained the electric field outside the shell (up to the coefficients c±), the 
next step is to match it to the inside solution at u = Ug using the matching conditions of 
eq. (2.5). To this end, we note that for u > Ug, the transverse EoM becomes 



El{u) + 



ijj 



fmf 



JmU 



E-iu] 



0, 



(3.10) 



where fm is the constant factor defined in eq. (2.4), originating from the relation between 
the inside and outside frequencies. The solution of this equation is easily found and can 
be expressed in terms of Bessel functions in the form 



E-{n) 



u\ J 



2Ju 



iY^ 



2Ju 



Jm 



(3.11) 



which we plug into a relation for the ratio C-/c-^ solved from the matching conditions, 



c_ 



-E'in duE. 



'TmE- duEii 



Eout duE- 



' fmE- ^rEo^^t 



(3.12) 



With this parameter determined, we have found the behavior of the (transverse) electric 
field for all values of u up to an irrelevant overall normalization constant. 



Repeating the above procedure for the longitudinal component of the electric field, we 
obtain the retarded Green's function of the photon upon an application of the relation 



R{uj, q, Us) = ^_ .^,^^ ^^ . (3.14) 



n(c2^,<7) = -^limj;^ (3.13) 

8 M->0 rj(u) 

separately to both the transverse and longitudinal components of Il^u{Q). Inserting the 
results to eq. (3.1), this gives us the trace of the corresponding spectral density Xm('^' ^' ""«) 
as well as its relative deviation from the equilibrium limit (denoted by 'th' below), 

XmI'^' g>^s) - (X th)M(^> 9) 
(Xth)/l(w, 

The behavior of this quantity — in particular the effect of the virtuality parameter c on 
its cj-dependence — is one of the most central issues to be studied in the following sections 
of this paper. 

3.2 Finite coupling corrections 

To study photon production outside the infinite 't Hooft coupling limit in principle ne- 
cessitates working out all the contractions in eq. (2.8) and subsequently determining the 
7-corrected gauge field EoM. This rather tedious exercise can, however, be avoided in the 
case of a transverse vector field, for which all the necessary formulae are available in the 
literature [28, 29, 40]. Taking advantage of these results, we will now study the behavior of 
retarded correlators of a transverse electric field E_\_{u) in the falling shell setup, keeping 
track of all terms linear in 7. It is, however, good to recall that according to eqs. (3.1)- 
(3.3), these results can be related to the photon production rate only in the limit of zero 
virtuality or c = 1. 

Following the treatment of [28], we begin by redifining the photon field according to 

^{u) = ny)E^{u), S(n)-i = l/v/7M + 7P(n), (3.15) 

where the function p{u) is given by 

, , ^2 (11700 -u^ [343897 -u(87539'u + 37760 ^2)]) 
viu) = — ^ , -^ . (3.16) 

288v7R 
This reformulation has the advantage of reducing the system to a simple Schrodinger-type 
problem with the action 






^^c^ + du^ 



(3.17) 

where the boundary term has the form ^(u) = ^'(n)^(n) and the EoM for ^(u) reads 

^"(u) - F(n)^(u) = 0, (3.18) 



1 

+1838319n'^ + ^^(4770 + 11700u2 - 953781^^ + 1011173n'') 



V(u) = --^ ^^ ^^ - 7^( - 11700U + 2098482n3 - 4752055^^ 

V ^ f2 „. '144\ 



-ui 



2(4770 + 28170'u2 - 1199223n*^; 



(3.19) 



From this point on, solving for the retarded Green's function of the photon proceeds 
along the very same lines as in the A = oo limit. Outside the shell, we write the ingoing 
and outgoing solutions to the EoM in the forms 

%n/out{u) = (1 - u^f (^L%(u) + l4^/.Jn) + 0(7^)) , (3.20) 

then separate eq. (3.18) into 0(7*^) and 0(7""^) parts, and finally solve "^^^'{u) and "^^^'{u) 
using the same boundary conditions as before. Just like at 7 = 0, this gives us the outside 
solution in terms of two unknown constants, 

*+(n) = C+^in{u) + C_^out(w) . (3.21) 

Inside the shell, we on the other hand use the fact that the metric (2.11) does not receive 
7-corrections [41, 42], implying that the 0(7) part of the ^_ field, 

^_ (n) = ^L°^ {u) + 7^L^) (n) , (3.22) 

originates solely from the relation (2.4) between the inside and outside frequencies, where 
we must simply change 

fm^ fl = fmK\us) . (3.23) 

This in particular means that ^_ (n) can be read off directly from eq. (10) of [23] (noting 
that inside the shell '^{u) = E±{u)), while ^_ (n) is available through a simple expansion 
of Bessel functions. 

Moving on to patching the two solutions together, we note that when written in terms 
of the field ^, the matching conditions of eq. (2.5) take the forms 

^.(us) = \/7I^+(n,)/S(n,) , (3.24) 

from which one can derive the 0(7) equivalent of eq. (3.12). This leads to a 7-correction to 
the ratio c_ /c+ , with which the transverse component of the spectral density is available 
from the relation [28] 



x±{^,q,us,v = -^^— 1 — ^7 Im 



expanded to linear order in 7. This result is finally used to define the quantity 



(3.25) 

M=0 



R±{u},q,u,,-f} = \ f-^ - \ ' (3-26 

ixth)±[^,q,7) 

which serves as our 0(7) generalization of eq. (3.14) and the behavior of which we will 
thoroughly analyze in the next section. There we will also compare the large-w limit of our 
numerical results to an analytic WKB expression derived for the c = 1 spectral density in 
appendix A. 



4 Results 

After explaining the main steps of our computation above, we will now move on to a nu- 
merical investigation of the corresponding results. This will be performed in three parts. 
We begin from the limit of thermal equilibrium, and inspect the effect of infinitesimal per- 
turbations on the system by solving the QNM spectrum from the thermal photon correlator 
and by analyzing the flow of the poles as functions of A. Next, we move to the falling shell 
setup, and perform a study of the virtuality dependence of the out-of-equilibrium spectral 
density of the photon, first in the A = cxd limit and then taking the leading 7-corrections 
into account. Finally, we specialize to the UV (a; — )• 00) limit of the 0(7) spectral density 
of real photons, and inspect how well the analytic WKB solution presented in appendix A 
approximates our full numerical results. 

4.1 Quasinormal mode spectrum 

Quasinormal modes (in equilibrium) are the strong coupling equivalent of the quasiparticles 
encountered in weakly coupled thermal field theory, and characterize the response of the 
system to infinitesimal external perturbations. On the dual gravity side, their spectrum is 
obtained by considering linearized fluctuations of some bulk field with infalling boundary 
conditions at the black hole horizon, and solving for the poles of the corresponding retarded 
Green's function. The QNM spectrum has the generic form 

ujniq) = Mn{q) - iTn{q), (4.1) 

where q denotes the three-momentum of the mode, and the real quantities Mn{q) and r„(g) 
correspond to the energy and width (decay rate) of the excitation. 

For the case of the AA = 4 photon correlator considered in this work, the QNM spectrum 
obtains a simple form in the respective limits of g — t- and A = 00 [27, 39], 

ojniq = 0) = 2TTTn{±l - i) , n = 0, 1,2, ... (4.2) 

The fact that the widths of the excitations grow linearly as functions of their energy 
highlights the contrast between this system and one composed of longlived quasiparticles 
with lin{ojn) ^ Re(a;„). This observation can additionally be connected to the top-down 
nature of holographic thermalization at infinite coupling: Fluctuation modes with the 
highest energies are also the most shortlived ones, indicating that the hardest scales are 
the first ones to reach thermal equilibrium. 

In light of the above discussion, it is clearly very interesting to investigate, how the 
weak coupling quasiparticle spectrum turns into the QNM one, as the system becomes 
more strongly coupled. With the formalism introduced in the previous two sections, this 
problem can be approached from the strong coupling side by determining the poles of the 
thermal photon Green's function including its leading finite coupling corrections. Recalling 
eq. (3.13), this calculation reduces to finding the zeros of the functions ip^iu = 0) and 
'i/'-jj (n = 0), introduced in eq. (3.20), which we write in terms of Probenius expansions [43] 

N N 

V^(°)(n,a;) = ^ a„(u;)(l - «)" , 4l\u,oo) = J^ 6„(a;)(l - n)" . (4.3) 

n=0 n=0 
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Figure 1. The flow of the four lowest q = i) quasinomial poles on the complex lo plane as functions 
of A. The first, second, third and fourth poles are denoted by circles, squares, diamonds and upward 
triangles, respectively, while the sole downward triangle corresponds to the fifth pole at A = 5000. 
The dashed lines of different colors correspond to 't Hooft coupling values ranging from A = oo to 
A = 1000, and have been drawn here merely to quide the eye. 
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Figure 2. The flow of the real and imaginary parts of the flrst poles of flg. 1 (circles and squares) 
as functions of A. 

The coefficients a„ and 6„ appearing here are determined via an analytic solution of recur- 
sion relations provided by the EoM (3.18), while the parameter A^ is chosen large enough so 
that the behavior of both functions becomes stable. The QNM spectrum is then obtained 
by making an ansatz for the frequency, cj„ = a;„ + '^ujn , substituting this to the equation 
il>ln {0, (jj) +7^^111 (0, oj) = 0, and finally solving the system numerically order by order in 7. 
The results of the computation just described are displayed in figs. 1-3. In the first of 
these, corresponding to the q = case, we show the flow of the first four QNM poles on 
the complex cj plane, of which the first two are further inspected in fig. 2. As the 't Hooft 
coupling is decreased from the A = 00 limit towards A = 1000, we observe a clear bending 
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Figure 3. Same as fig. 1, but for g = 1, i.e. q = 27rT. 

of the spectrum away from the Unear form. In particular, we witness a rapid increase of 
the imaginary parts of a;„, as A is lowered, as well as a strong dependence of the magnitude 
of the 0(7) term on the index n. These observations are suggestive of the QNM spectrum 
moving towards a quasiparticle one at smaller couplings, and are clearly consistent with 
the weakening of the top-down pattern of the thermalization process, cf. [30]. It should be 
noted, though, that our strong coupling expansion can be trusted to converge only when 
the relative deviation of a pole location from its A = 00 value is small, which clearly is not 
the case for all of the poles shown in fig. 1. Finally, in fig. 3 the same exercise is performed 
for the case of q = 1 (or q = 2ttT) with qualitatively similar results. 

4.2 Virtuality dependence of the spectral density 

Next, we move on to study the behavior of the retarded photon correlator in a specific 
out-of-equilibrium setting, namely during the thermalization process corresponding to the 
collapsing shell model. Beginning from the limit of A = 00, we want to complement 
our earlier work on real and maximally virtual photons [22, 23] by letting the virtuality 
V = [up' — c^^juP' be a free parameter. To this end, we parameterize q = cu and plot in fig. 4 
the behavior of the spectral density x^ and its relative deviation from the thermal limit, R 
(cf. eq. (3.14)), for c = 0, 0.8 and 1. In analogy with the observations of [22, 23], we again 
witness an oscillation of the off-equilibrium spectral densities around the corresponding 
equilibrium values, with the amplitude of the oscillations decreasing as the shell approaches 
the horizon radius.^ The most important effect of the virtuality is clearly visible in the 
behavior of the relative deviation R: The larger the virtuality of the photon is (the smaller 
the value of c), the smaller is the amplitude of the oscillations. This implies that maximally 



^ Since the other effects of the shell location are very minor, we have in ah of our plots set this parameter 
to the rather arbitrary value Ms = 1/1.1'^. 
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Figure 4. Left: The A = oo spectral density for photons of different virtuahties in thermal 
equilibrium (dotted lines) and in an off-equilibrium state described by the collapsing shell model 
with Us = 1/1.1^. Right: The corresponding relative deviations from the thermal limit, R, for c — 
(solid blue line), c = 0.8 (dotted red line) and c = 1 (dashed black line). 
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Figure 5. As in fig. 4, but for A = 300 (left) and A = 100 (right). 

virtual photons, i.e. dileptons at rest, are the first ones to thermalize, confirming the 
observations made in somevi^hat different contexts in [44, 45]. Finally, we note that in all 
cases studied the amplitude of the oscillations in R decreases at large w, consistent with 
the known top-down nature of thermalization at infinite 't Hooft coupling. 

At finite values of A, one expects to witness a qualitative change in the behavior of the 
deviation function (now R±), as reported for real photons (c = 1) in [30]. To investigate 
what happens at nonzero virtuality, we determined the transverse spectral density x± and 
the corresponding R± for the same values of c considered above, but setting now A = 300 
and 100, cf. fig. 5. Similarly to the case of real photons, we again observe the asymptotic 
behavior of the fluctuation amplitude changing from a 1/w suppression towards linear 
growth as the coupling is decreased. This behavior, however, appears to depend on the 
virtuality rather strongly, with the transition happening later as c is decreased; in the case 
of maximal virtuality, c = 0, the amplitude of R± even appears to approach a constant at 
large u. Finally, we note that the A = oo observation of the amplitude of R± decreasing 
with increasing virtuality seems to hold at all coupling strengths considered. 

Just as in the case of [30], one must exercise some caution with the above results, 
not least because both the quasistatic approximation and the strong coupling expansion 
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Figure 6. A comparison of our exact numerical result for the relative deviation R (dashed black 
line) with its analytic asymptotic limit derived from eq. (4.4) (solid red line), with both results 
evaluated using Ug — 1/1.1^ and 7 = 100. 



we have employed have finite regions of applicability. We note, however, that in the pa- 
rameter ranges considered above, our results should be soHd. First of all, a comparison of 
the relevant frequencies to the inverse time scale associated with the motion of the shell, 
cf. eq. (2.1), clearly indicates that neglecting the latter is a reasonable approximation to 
make. And second, it can be explicitly checked that with the exception of the very largest 
frequencies studied, the deviation of x± from its A = 00 limit is always at most of the order 
of 10%, giving us confidence that the strong coupUng expansion is indeed well behaved. 

4.3 UV limit of the spectral density 

As the final test of our results, we want to study the extent one can understand the 
UV (large w) behavior of the spectral density of real photons analytically, recalling of 
course that in this formal limit the strong coupling expansion most likely breaks down. 
The calculation is nevertheless highly relevant for the arguments of [30] concerning the 
asymptotic behavior of R being indicative of the thermalization pattern of the system. 

In appendix A, it is shown that a WKB-type calculation for the retarded photon 
correlator results in a relatively compact result for the large-w behavior of the 0(7) spectral 
density. 
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Here, one should note that while the order 7*^ off-equilibrium terms residing on the second 
line behave at large frequencies like 1/w, the corresponding 0(7) terms on the third line 
indeed exhibit a linear behavior in lo, consistent with our observations. In fig. 6, this 
comparison is made more quantitative by displaying the relative deviation R obtained from 
this result with our full numerical calculation. As can be seen from the figure, the agreement 
between our numerics and the asymptotic large-o; limit is quite impressive already at 
moderate frequencies, improving consistently as one moves to even higher values of w. 

5 Discussion and conclusions 

5.1 Summary of results 

In the paper at hand, we have addressed a number of open issues related to the production 
of (prompt) photons in a thermalizing holographic plasma, motivated by our earlier works 
on the same subject [22, 23, 30]. Our initial studies in the infinite 't Hooft coupling 
limit [22, 23] only covered the two limiting cases of real and maximally virtual photons 
(c = 1 and c = 0, respectively), and thus left open the determination of the full virtuality 
dependence of the results. More importantly, [30] made some highly nontrivial claims 
about the thermalization pattern of the holographic plasma shifting away from the usual 
top-down behavior, as the 't Hooft coupling is decreased from the A = 00 limit. This 
conclusion was mostly based on the UV behavior of the retarded photon Green's function, 
which was, however, studied only numerically. Moreover, all of these results were derived 
using only one particular model of holographic thermalization, and a clear need for a more 
universal analysis was identified. 

In the following, we briefly summarize the findings we have made in the present work, 
and explain how we believe they shed light on the issues identified above: 

Quasinormal mode analysis 

In section 4.1, we investigated the flow of the photon QNM spectrum, as the first strong 
coupling corrections were added to the A = 00 limit a;„ = n(ibl — i), and the 't Hooft 
coupling was decreased. Our findings, displayed in figs. 1-3, are strongly suggestive of the 
bending of the spectrum upwards on the complex uj plane. One must, however, naturally 
exercise some caution with the exact locations of the poles that have shifted an 0{1) 
distance from their A = 00 values, as such a shift raises doubts on the convergence of the 
strong coupling expansion employed. 

If taken at face value, our QNM result is significant for two distinct reasons. First, it 
indicates a clear shift of the system towards one described by longlived quasiparticles, often 
encountered in weakly coupled thermal plasmas. The fact that the widths of the excitations 
are no longer proportional to their energies can also be interpreted as a sign of the weakening 
of the top-down thermalization pattern of the infinitely strongly coupled system, and thus 
strengthens the arguments made in [30]. Second, perhaps more importantly our QNM 
study is completely independent of the thermalization model used in our earlier works, 
and any approximations associated with it. This result is thus one of considerably more 
universal value than those relying on the collapsing shell setup. 
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Virtuality dependence 

In section 4.2, we analyzed the virtuality dependence of the photon spectral density both 
in the strict A = oo limit and including the leading 0(7) corrections. At infinite coupling, 
the results depicted in fig. 4 confirmed the observations of [44, 45], suggesting that highly 
virtual field modes are the first ones to thermalize. At the same time, fig. 5 extended these 
results to finite coupling, revealing a strong virtuality dependence in the large-w behavior 
of the 0(7) spectral density. In particular, with maximal virtuality the relative deviation 
of the spectral density from its thermal limit was seen to approach a constant in the a; — )■ cxd 
limit, in contrast with a linear rise observed for all nonzero values of c. 

UV behavior of the spectral density 

Section 4.3, as well as appendix A, finally dealt with the analytic evaluation of the large-w 
limit of the photon spectral density in the limit of zero virtuality but including strong 
coupling corrections. A WKB-type calculation performed in the appendix was seen to 
result in a rather compact expression for the formal UV limit, displayed in eq. (4.4), 
which furthermore agreed with our full numerical result already at intermediate frequencies, 
cf. fig. 6. This result clearly strengthens the conclusions made in [30], as it confirms the 
linear increase of R in the w — t- 00 limit. At the same time, it is nevertheless good to recall 
that due to the leading 7-corrections to R being accompanied by a factor of lo, this formal 
limit will never be reached within the validity of the strong coupling expansion. 

5.2 Future directions 

The most important qualitative outcome of the present paper is clearly the strengthening 
and confirmation of the most central conclusions of [30]. Both the generalization of our 
earlier calculations to nonzero virtuality and the analytic treatment of the large-w limit of 
the photon spectral density point towards the weakening of the top-down thermalization 
pattern as the 't Hooft coupling is decreased from the A = 00 limit. At the same time, 
the determination of the photon QNM spectrum revealed compelling evidence for the fiow 
of the system towards one described by quasiparticles, as the coupling is decreased. As 
discussed above, this pattern is furthermore in accordance with the weakening of the top- 
down thermalization behavior, as it points towards the relaxation of the correspondence 
between the energy and inverse lifetime of the field fiuctuations once the value of the 
coupling is lowered. 

While the results obtain in this paper are in many ways promising, much more remains 
to be done. On one hand, progress in the study of fully dynamical thermalization away from 
the infinite coupling limit would clearly be extremely valuable, generalizing results such as 
[5, 11, 45-50] to finite A. At the same time, our current results on photon production should 
be straightforward to generalize e.g. to the case of energy momentum tensor correlators, 
providing direct information on the thermalizing plasma constituents. Both of these topics 
are indeed amongst the issues we are currently working on. 
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A WKB determination of the spectral density 

In this appendix, we will perform an analytic WKB-type study of the large-w behavior of 
the retarded Green's function of a real (transverse) photon field, keeping all terms to linear 
order in 7. This exercise is greatly simplified by the fact that the EoM (3.18) satisfies the 
conditions of Theorem 3.1 of [51], concerning the asymptotic behavior of the solutions of 
differential equations of a certain type. Applied to the computation at hand, this implies 
that just like in the 7 = case studied in [23], our general large-w solution can be expressed 
in terms of a linear combination of Airy functions, 

^+(u) = h-'/\-u)^ci AiiCo^/'^Ci-u)) + C2 Bi((2;2/3^(_^))| . (a.i) 

Here, the q denote unknown coefficients (with no relation to the previous c±), while the 
functions h{—u) and C,{—u) read 




1^ 

g{u) = 2arctanh(\/u) — 2 arctan(-v/n) 

^3/2 / \ 

-7 422730 - 2699862n2 - 7078211n^ , 

' 5544 V ) 

u { I- V? f 

hi-u) = --, ^rrr- -I 1 + 7 — -— 16470 - 245442tt2 - 1011173u^ 

(1 -v?YC{-u) 1 144 V 

Similarly to the corresponding 7 = expressions of [23], the function g{u) is real and 
positive in the range < u < 1 for all reasonable values of 7 (this ceases to be true for 
A < 4.2, but such small values of the coupling are not considered in this work). 

Following the steps explained in [23], it is straightforward to verify that the (correctly 
normalized) linear combinations of the Airy functions satisfying infalling and outgoing 
boundary conditions at the horizon take the respective forms 

^in W = ^^-^^P^ Ai{u'/\{-u)) , (A.5) 

^out(n) = ^K"}-""^ UKi{u:^"a-u))+Bi{u^/\{-u)) 
2aOn„t V 



'out 
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where we have denoted 

2(i/2g-i7r/4 



a = ; ;:. , /? = 2-e---/^ , (A.6) 

An = 1 - TT^ (9419333 + 9355343z^ + 1453760cuM , (A.7) 

11088 y J 

/3out = 1 ^ (9419333 - 9355343^^ + 1453760wM . (A.8) 
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Using the identities 
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,^^^,+0{y-'/^), (A.9) 

9„-2J/3y3/2 -^i/S 2/3/2 

^-Wi2Bi(e-/3,) ^ ^^ ^ + 0(y-^/^) (A.IO) 

for the asymptotic behavior of the Airy functions, it is then rather straightforward to derive 
the large-cj limit of the retarded Green's function, 

/V27^2ol/3T^,'2/Q^ i_22i<ig-i7r/2('i+l/3) ft2_£^ 
IlT'^{LO,Us,l) = - ' ,^f'^\ l + b^) A^^li:±(_^)2/3 (A.ll) 

-L '^ ' '^' '^ 8 r 1/3 ^ '^1, 22iie-^'^/2('^-i/3)-fii^£z^^ ' '^ ^ 

Pout C+ 

where the ratio of c_ and c+ remains as the only unknown parameter. 

In determining the large-w limit of C-/cj^^ we finally take advantage of eq. (3.12), to 
which we plug the above asymptotic forms of the infalling and outgoing field modes. After 
some algebra, this produces the result 

C./C+ = Co + 7Ci + 0(72), (A. 12) 

Co „^ I with2-2^'^n;3/2^-i/r^^(l + 7r^^)e*'^("/2+9oK))^ (A.13) 



w— >oo 



Ci _ ^ ^^ 2-2-nI/2(l _ ^2)^g^<i^/2+3oK)) ^ (^.U) 
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where go(^) = 5(^)1^=0. Inserting these expressions to eq. (A. 11), we then obtain an 
analytic, albeit quite complicated expression for the retarded correlator. Upon taking the 
imaginary part of this result and dropping all terms beyond linear order in 7, this function 
however simplifies considerably, producing the rather elegant 0(7) generalization of eq. (35) 
of [23], 



iV2r2 35AT(2/3).2/3 
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Oh')]. (A.15) 
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In this expression, the first Hne represents the equiUbrium hmit of the result, while the 
second contains the 0(7'') and the third the 0(7^) non-equilibrium corrections to it. 
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